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Abstract 

This paper is devoted to study the charged perfect fluid cylindrical 
gravitational collapse. For this purpose, we find a new analytical solu- 
tion of the field equations for non-static cylindrically symmetric space- 
time. We discuss physical properties of the solution which predict 
gravitational collapse. It is concluded that in the presence of electro- 
magnetic field the outgoing gravitational waves are absent. Further, 
it turns out that when longitudinal length reduces to zero due to re- 
sultant action of gravity and electromagnetic field, then the end state 
of the gravitational collapse is a conical singularity. We also explore 
the smooth matching of the collapsing cylindrical solution to a static 
cylindrically symmetric solution. In this matching, we take a special 
choice of constant radius of the boundary surface. We conclude that 
the gravitational and Coulomb forces of the system balance each other. 
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1 Introduction 



Gravitational collapse of a massive star occurs when all the thermonuclear 
reactions in the interior of a star could not favor the pressure against grav- 
ity. Gravitational collapse is one of the most important problems in general 
relativity. The singularity theorems ^ state that there exist spacetime sin- 
gularities in the realistic gravitational collapse. To investigate the nature of 
spacetime singularity, Penrose 2 > suggested a hypothesis known as Cosmic 
Censorship Hypothesis (CCH). It states that the final fate of gravitational 
collapse of a massive astrophysical object is always a black hole. This is 
equivalent to saying that the singularities appearing in gravitational collapse 
are always clothed by an event horizon. 

Many attempts predicted that final fate of gravitational collapse of the 
massive star might be a black hole or naked singularity depending upon 
the choice of initial data. In this chain, Virbhadra et al. 3 ) introduced a 
new theoretical tool using the gravitational lensing phenomena. In a recent 
paper A \ Virbhadra used the gravitational lensing phenomena to find an 
improved form of the CCH. The classical paper of Oppenheimer and Snyder 
5 ) is devoted to study dust collapse according to which singularity is neither 
locally or globally naked. In other words, the final fate of the dust collapse is 
a black hole. Many people 6 ) extended this work for physically existing form 
of fluid with cosmological constant in spherically symmetric background. 

In order to generalize the geometry of the star, people worked on grav- 
itational collapse using cylindrical symmetry. The existence of cylindrical 
gravitational waves provides a strong motivation in this regard. Bronnikov 
and Kovalchuk 7 ) were the poineers to the work on gravitational collapse with 
non-spherical symmetry. Later on, the same authors 8 ) extended it for some 
non-spherical exact model. Nolan 9 ) investigated the naked singularities in 
the cylinderical gravitational collapse of counter rotating dust shell. 

Hayward 10 ) studied gravitational waves, black holes and cosmic strings in 
cylindrical symmetry. Sharif and Ahmad n ) analyzed cylindrically symmetric 
gravitational collapse of two perfect fluids using the high speed approxima- 
tion scheme. They investigated the emission of gravitational radiation from 
cylindrically symmetric gravitational collapse. Di Prisco et al. discussed 
the shear free cylindrical gravitational collapse using junction conditions. 
Nakao et al. 13 ^ studied gravitational collapse of a hollow cylinder composed 
of dust. 

Gravity is the weakest interaction among all the natural forces. The be- 
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havior of electromagnetic field in gravitational field has been the subject of 
interest for many years. Thorne 14 -* developed the concept of cylindrical en- 
ergy and investigated that a strong magnetic field along the symmetry axis 
may halt the cylindrical collapse of a finite cylinder before it reached to sin- 
gularity. Oron 15 ) studied the relativistic magnetized star with the poloidal 
and toroidal fields. Thirukkanesh and Maharaj found that the inclusion 
of electromagnetic field in gravitational collapse would counterbalance the 
gravitational attraction by the Coulomb repulsive force along with pressure 
gradient. In recent papers 17 )- 19 ), we have investigated the effects of elec- 
tromagnetic field on the perfect fluid collapse by using junction conditions 
in spherically symmetric background with positive cosmo logical constant. 
It has been found that electromagnetic field reduces the pressure and favors 
the naked singularity formation but cannot play a dominant role. Thus black 
hole was formed as a final state of the gravitational collapse. 

In this paper, we study the cylindrically symmetric charged perfect fluid 
gravitational collapse. The main objective of this work is to study the final 
fate of charged perfect fluid gravitational collapse in the cylindrically sym- 
metric background. The plan of the paper is as follows: In the next section, 
we discuss the solution of the Einstein-Maxwell field equations. The physi- 
cal properties of the solution are discussed in section 3. Section 4 gives the 
derivation of the matching conditions. We summarize the results in the last 
section. 

Geometrized units (i.e., the gravitational constant G—l and speed of light 
in vacuum c = 1) are used. All the Latin and Greek indices run from to 3, 
otherwise, it will be mentioned. 

2 Solution of the Einstein Field Equations 

This section is devoted to the solution of the Einstein field equations coupled 
with the charged perfect fluid as the source of gravitation distributed per 
unit length of the cylinder. The general cylindrically symmetric spacetime is 
given by the following line element 12 -* 

ds 2 _ = A 2 {dt 2 - dr 2 ) - B 2 d9 2 - C 2 dz 2 , (2.1) 

where A, B and C are functions of t and r. Here we take the following 
restrictions on the coordinates in order to preserve the cylindrical symmetry 
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of the spacetime 

-oo<t<oo, r>0, 0<9<2n, -oo < z < oo. (2.2) 

The proper unit length of the cylinder for the line element ( 12. ip is defined by 

I = 2nBC. (2.3) 
The Einstein field equations are given by 

G; = k(T^ + T^). (2.4) 
The energy-momentum tensor for perfect fluid is 

t; = ( P + P )ux - ps;, (2.5) 

where p is the energy density, p is the pressure and = A5® is the four- 
vector velocity in co-moving coordinates. The energy-momentum tensor for 
the electromagnetic field is given by 

r; (em) = ^(-F uX F»x + -fiF^F^), (2.6) 

where is the Maxwell field tensor. Now we solve the Maxwell's field 
equations 

Fpv = fa* - <f>^, F^, u = AnJ^ (2.7) 

where M is the four potential and J M is the four current. In co- moving 
coordinate system, the charge per unit length of the cylinder is assumed to 
be at rest so that the magnetic field will be zero. Thus we can choose the 
four potential and four current as follows 

</y =(</>(*, r), 0,0,0), J» = au^ (2.8) 

where a is charge density. The only non-zero component of the field tensor 
is 

F i = -F 10 = -|t (2.9) 
Thus the Maxwell field equations take the following form 



where dot and prime indicate derivatives with respect to time t and radial 
coordinate r, respectively. Integration of Eq. fl2.10p implies that 



where q(r) = 2n J a(ABC)dr, being the consequence of conservation law of 
charge, i.e., Jft = is known as the total amount of charge per unit length 
of the cylinder. It is mentioned here that Eq.( l2.11j) is identically satisfied by 
the solution of Eq. ( I2.10p . We can write Eq. 02. 9p as follows 

Fox = ~F 10 = (2.13) 
The non-zero components of T^ em ^ turn out to be 

rpO( em ) _ rpl(em) _ _rp2( em ) _ _rp3( em ) _ ^ ^ 



2tt(BC) 2 ' 
The electric field intensity is defined by 

£ < r - ( > = MW (2 ' 14) 

We assume that the charged perfect fluid distributed per unit length of the 
cylinder follows along the geodesies in the interior of the cylindrical symme- 
try. This requires that velocity should be uniform and acceleration must be 
zero which is only possible if A is constant and in particular, we take A = 1 
(for simplicity). Thus the field equations (12 .4p take the following form 
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We note that there are five equations and five unknowns B, C, p, p and E, 
thus we can find a unique solution. 

For this purpose, we adopt the method of separation of variables. The 
comparison of Eqs.f l2.18j) and (12. 19 j) give 

C G" B B" , nnn . 
-C + ^ = ~B + ^ (2 ' 20) 

which yields the necessary condition for pressure to be isotropic. We take 

B = f(r)g(t), C = h(r)k(t). (2.21) 

Using Eq.(j2j2H) in fT2TT6|) . we get 

f = ah L , k = 5g~ L , (2.22) 

where L 0, for non-trivial solution) is a separation constant while a and 
S are integration constants. Using Eq. (12.221) in (12.201) . it follows that 

9 k f" h" 

Since both sides are functionally independent, we put them equal to constant 
say M(^ 0) 

g k f" h" 

- - T = M = J —- — . 2.24 
g k f h 

Application of Eq.J222J to Eq. fl2~24l) leads to 

g g 2 _ M h" h' 2 M 

The solution to these equations is 

g (t) = f3 cos^L{Wt + t ), h(r) = p lC oshTTL(Sr + r ), (2.26) 

where j3o, to and r are constants of integration. Further, W and S are 
given by the following relations 



Using Eq.( g3SD in (127221) . it follows that 

= /3 2 cos^(H^ + t ), /(r) = & cosh ttt (Sr + r ), (2.28) 

where and are constants of integration. Thus the metric coefficients, 
given by Eq. (12.211) . turn out to be 

B = Vlcosh^ (Sr + r ) cos^ (Wt + to), (2.29) 
C = ^coshTTr (Sr + r ) cos^ (Wt + t„) , (2.30) 

where Q = /3o03, $ = fiifa- Consequently, the spacetime (2.1) takes the 
form 

ds 2 _ = dt 2 - dr 2 - ft 2 cosh^ (Sr + r ) cos^ (Wt + t )d6 2 

- ^cosh^rr (Sr + r ) cos"^ (Wt + t )dz 2 . (2.31) 

Using the following transformations 

S r f = Sr + r , Wt' = Wt + t , 9' = Q6, z' = Vz, 

this metric reduces to 

ds 2 _ = dt' 2 -dr' 2 -cosh^(Sr')cos^(Wt')d6' 2 

- coshTTL(Sr')cos^(Wt')dz' 2 . (2.32) 

By assuming Q — 1, it is clear that the above metric preserves cylindrical 
symmetry with the restriction on coordinates given by Eq. fl2.2p . Here we 
take 

B = cosh ttt (Sr') cos^r (wt') , C = cosh^ (Sr') cos^r (Wt') . 

3 Physical Properties of the Solution 

Here, we discuss some physical and geometrical properties of the solution. 
The physical parameters, i.e., pressure p, density p, and the electric field 
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intensity E for the metric (I2.32p are given by 



1 S 2 Atan 2 (Wt')W 2 L W 2 
P ~ 1&7T (1 + L) (1 — L) 2 + (L-1)(2L- 

1 2L{l + L) 2 W\ec 2 {Wt')-{l + LfW 2 -{L-lfS 2 
L 32^ (1 — L 2 ) 2 

2LSHec 2 h{Sr>) , 
+ (L + l) 2 IJ ' ^ 
1 r -g 2 (l + L + L 2 + Lsec 2 h{Sr')) LW 2 tan 2 {Wt') } 

9 ~ (lTT? + (L-iy 1 {3J} 

We would like to mention here that Eqs.([2]29D, ( ESQ} and (EHJ) - (J3l3]) satisfy 
all the field equations with the restriction on constants given by Eq. (12.27]) . 
The proper unit length of the cylinder for the new metric is given by 

I = 2nBC = 27rcosh(,Sr , )cos(iyt') (3.4) 
and the longitudinal length in this case is 

l = BC = cosh(Sr')cos(Wt') . (3.5) 
The rate of change of longitudinal length is 

1 = -Wcosh(Sr')sm(Wt'), (3.6) 

where negative sign shows that motion is directed inward. 20 > Thus such mo- 
tion represents gravitational collapse of the charged perfect fluid distributed 
per unit length of the cylinder. 

In order to analyze the nature of singularity of the solution, we use the 
curvature invariants. Many scalars can be constructed from the Riemann 
tensor but symmetry assumption can be used to find only a finite number of 
independent scalars. Some of these are 

r> D „ab r> td D paft p p pafecd r> pa6 rycd 

il\ — it — g K a b, rt2 — K a bK , H3 — K a bcd-tl , -K4 — J^cd K ab- 

Here, we give the analysis for the first invariant commonly known as the 
Ricci scalar. For the metric (I2.32p . it is given by 

R = l(BC -b6- B"C - C"B + BC- B'C'), (3.7) 
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where / is given by Eq. fl3.5p . We see that Ricci scalar as well as all the other 
curvature invariants and physical parameters of the solution are finite for 
r' — > 0. Thus r' = is the conical singularity of the metric f)2.32p . 

Now we analyze the values of the constants for which the solution is 
physical. In this solution, L and M are non-zero separation constants for the 
non-trivial solution while the rest are integration constants that are removed 
by applying the transformations to Eq. fl2.3ip and by evaluating the physical 
parameters from the field equations. From Eq.f l2.27j) . it is clear that the 
constants W and S cannot be chosen arbitrarily. These are non-zero because 
M 7^ for non-trivial solution. Further, for W and S to be real, there are 
following four possible solutions: 

1. L<-1, M>0; 2. L>-1, M < 0; 
3. L>1, M>0; 4. L < 1, M < 0. 

Keeping in mind these restrictions on the constants, we find that the cases 
1 and 2 lead to non-physical solutions (i.e., negative energy density for the 
arbitrary choice of coordinates). In the case 3, for < M < 0.5 and 1 < 
L < 1.9, there exists a physical solution which represents the gravitational 
collapse. The graphs 1-4 in this case indicate that all the physical quantities 
become homogeneous. Thus the geodesic model with charged perfect fluid 
distributed per unit length of the cylinder is free of initial inhomogeneities. 

It is interesting to mention here that pressure remains function of time 
only for the geodesic model that is analogous to the spherical case. 21 ) In the 
case 4 for 0.63 < L < 0.95 and —1 < M < —0.10, all quantities except 
pressure behave like the case 3, in this case pressure is negative indicating 
a dark energy solution. As long as the realistic energy condition p + 3p > 
holds, the gravity remains attractive. 22 -* However, the violation of this 
condition i.e., p + 3p < due to negative pressure, leads to the repulsive 
gravitational effects. Thus in the relativistic physics, negative pressure acting 
as a repulsive gravity plays the role of preventing the gravitational collapse. 
We are interested to study the gravitational collapse which is the consequence 
of attractive gravity, so the case 4 is not interesting here. Thus the only 
interesting case is the case 3. 

Now we proceed to discuss the energy conditions for case 3 which are 
given by the following relations ^ 

1. Weak energy condition: 

P>0, p + p>0, (3.8) 
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2. Dominant energy condition: 

P + P>0, p-p>0, (3.9) 

3. Strong energy condition: 

p + 3p>0. (3.10) 

For the pressure and energy density given by Eqs. (13. 11) and (13.31) respectively, 
it follows that 

1 (l + L + L 2 + Lsec 2 M^)) 

P P 16iv [ l + L X 1 + L 1 

21**™)], (3.11) 

L-1 V 2L-1 L-l ;J ' V ; 

= i^ [ TTl {1 + 2 iTl } 

^. (M t^_ i )L (312) 



L-V L-l 2L-V 

p+3p _^_ 2 JLh£±^C*3) } 

W 2 , 1 tan 2 W^' 

10L— — )]. (3.13) 



L-1 V 2L-1 L-l 

Notice that all these equations satisfy Eqs.flBTHD- dBTTUD for < M < 0.5, 1 < 
L < 1.9, < t' < 1 and < r'. 

The rate of change of longitudinal length in Figure 1 shows that the 
longitudinal length is a decreasing function of time, thus the resulting solution 
represents the gravitational collapse. 23 - ) The collapse starts at some finite time 
and ends at t' — 1, where longitudinal length of the cylinder reduces to zero. 
Further, energy density is an increasing function of time shown in Figure 2. 
This is the strong argument for a model to collapse. The pressure in the 
interior of cylinder starts decreasing as shown in Figure 3. This causes to 
initiate the gravitational collapse, more matter is concentrated in the small 
volume, hence density goes on increasing. 

It is to be noted that decrease in the proper unit length of the cylinder, 
increases the interaction between the electric charges and a strong electro- 
magnetic tension inside the cylinder is created. This is an increasing function 
of time as shown in Figure 4. The coupled action of electromagnetic and grav- 
itational forces play a dominant role to reduce longitudinal length of cylinder 
to zero. 
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Velocity Graph 




Figure 1: Decrease in longitudinal length with the passage of time for < 
M < 0.5 and 1< L < 1.9 (Color online). 



Matter Density Graph 




Figure 2: Increase of density with the passage of time for < M < 0.5 and 
1 < L < 1.9 (Color online). 

P Pressure-Time Graph 
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Figure 3: Decrease in pressure with the passage of time for < M < 0.5 and 
1 < L < 1.9. 
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Electric Field Intensity Graph 




Figure 4: Increase of electric intensity with the passage of time for < M < 
0.5 and 1 < L < 1.9 (Color online). 

The nature of the collapse can be seen as follows: When latitudinal and 
vertical lengths of the cylinder reduce to zero, there is a complete collapse. 
From the metric (2.32), we have gee = B 2 , g zz = C 2 . Since singularity anal- 
ysis implies that the Ricci scalar diverges at a point where the longitudinal 
lengths I = BC = 0. Thus when the longitudinal length as well as the lat- 
itudinal and vertical lengths reduce to zero, we obtain a conical singularity 
at r' = 0. 

The conical singularity has effects on the gravitational collapse of non-zero 
mass objects. The conical singularity belongs to the class of line singularities 
which are gravitationally weaker than point singularities and stronger than 
the plane singularities. The tidal forces for point and line singularities are so 
strong that these can crush an object of finite size and non-zero mass to zero 
volume. However, it was pointed out by Nakao et al. 24 -* that only conical 
singularity is the exceptional line singularity which does not crush an object 
to zero volume that collapses onto it. Hence, it is concluded that massive 
objects of finite size are collapsed on it without crushing to zero volume. The 
reason of non-crushing to zero volume does not imply that tidal forces are 
weak but it may be due to the geometric structure of the conical singularity. 

4 Matching Conditions 

Following, 21, 25 > we proceed to cut the spacetime and match it with another 
spacetime, which represents the exterior region of the collapsing cylinder. We 
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match the charged perfect fluid solution with the electro-vacuum solution. 
For this purpose, we consider the Darmois junction conditions, 26 - 1 which 
require that the first and second fundamental forms (that are the line ele- 
ments and the extrinsic curvature respectively) must be continuous over the 
boundary surface S. 

We assume that the 3D timelike boundary surface £ splits the two AD 
cylindrically symmetric spacetimes V + and V~. The metric which describes 
the internal region V~ is the charged perfect fluid solution given by Eq. (12.32]) . 
For the representation of the exterior region V + , a charged cylindrically 
symmetric electro- vacuum solution is taken as 27 > 

ds\ = HdT 2 - —dR 2 - R 2 {d9 2 + dz 2 ), (4.1) 
H 

where H(R) = ^ — M and Q are the mass and charge per unit length 
of the cylinder, respectively. This choice of the exterior solution in V + region 
is compatible with the charged perfect fluid solution in the interior region 
V~ for their smooth matching over the boundary surface S. 

Now the boundary surface E in terms of interior and exterior coordinates 
can be described by the following equations 

/_(r',t')= r'-4 = 0, (4.2) 
U(R,T)= R-R s (T) = 0, (4.3) 

where r's is a constant. Using these equations, the interior and exterior 
metrics on E take the following form 

(dsl ) s = dt 12 - B 2 d6' 2 - C 2 dz' 2 , (4.4) 

(ds 2 + h = [H(Rx) - Jj^i^yT 2 - Rl(d9 2 + dz 2 ). (4.5) 

We assume g 00 > in Eq. (14.51) so that T is a timelike coordinate. 
The continuity of the first fundamental form gives 

(Bh = (Ch = R*, (4-6) 

1 dR 

[H(Rz) - -W{RV) { ^ )2]l2dT = {dt% - (47) 

The components of extrinsic curvature Kf- in terms of interior and exterior 
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coordinates are 



K, 



oo 



0, # 2 " 2 = #33 = (BBh, 




(^t T tt _ T t^tt T tt3 + 

2 ai? 



3 d# 
2HlR 



Tti?t 2 ) 




#3+ = (#i2Tt) E , 



(4.8) 



where dagger f and bar ~~ represent differentiation with respect to the new 
coordinates t' and r' respectively. The continuity of the extrinsic curvature 
components with Eqs. fl4.6p and (|4.7p leads to 



Here Eq. ( 14.9}) implies that the boundary surface £ represents a cylinder 
with constant proper unit length which behaves as boundary of the inte- 
rior charged perfect fluid distributed per unit length of the cylinder. Thus 
it connects the interior charged perfect fluid solution to the exterior electro- 
vacuum solution. Using this equation, Eq.f l4.10p reduces to M = (fg)s- The 
parametric representation of this equation implies that the gravitational and 
Coulomb forces of the system balance each other on the boundary surface E. 
This consequence in the absence of pressure on the boundary, is equivalent 
to the result found by Thirukkanesh and Maharaj. 16 -* 

5 Discussion 

In this paper, we find an analytical solution to the Einstein field equations 
coupled with the charged perfect fluid distributed per unit length of the cylin- 
der. It has been found that all the physical parameters become homogenous 
but isotropic pressure is function of time only. This property (i.e., pressure 
being only time dependent) of cylindrically symmetric geodesic model is sim- 
ilar to the spherical case. Further, all the energy conditions are satisfied for 
a range of separation parameters and coordinates for which the solution is 
physical. All sorts of singularities of the solutions are discussed in detail. 
It is found that physical singularity of the solution occurs at a point where 
proper areal radius of the cylinder reduces to zero. 



(5 f ) S = 0, 



(4.9) 



M=[X + ^( J Bt 2 - J B 2 )] E . 
2B 4 V 



(4.10) 
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The physical and geometrical properties of the solution such as increase 
of density and decrease in proper unit length of the cylinder with respect 
to time represents the gravitational collapse. Also, the time interval for the 
collapse has been investigated. It is found that the electric field intensity of 
the system increases with time. This implies that as longitudinal length of 
the cylinder decreases, charges come close to each other and the interaction 
between the charges is increased. 

In general, it is known 19 )> 28 ) that the presence of an electromagnetic field 
in a geometry causes to disturb its generic properties which may result in the 
form of oscillation of spacetime. Here, we discuss the absence of oscillation 
in geometry by showing the absence of outgoing gravitational waves. 

Following Pereira and Wang 29 \ the component of the Weyl tensor for our 
solution, is \t> = -C^ u x a L^M u L X M U = 0, where LP and M v are null vectors. 
It means that there does not exist outgoing gravitational waves implying that 
oscillations are absent in the geometry of the spacetime. Thus there is no 
energy loss and hence no bouncing. Further, the absence of fluctuations in 
the energy density graph 2 indicates the absence of bounce, oscillation and 
gravitational waves. The prediction that the fluctuation in energy density 
represents gravitational waves is recently made by Hussain et al. 30 ^. 

We would like to remark here that in this case electromagnetic field is 
weak as compared to matter field which is obvious from Figures 2 and 4. This 
condition [E 2 < p) for the weak electromagnetic field was stated by Tasagas 
2S ) and used by us to preserve the generic properties of the Friedmann 
universe models. In a paper 29 \ $o 7^ implies that there exist outgoing 
gravitational waves from the cylindrical symmetric radiating fluid gravita- 
tional collapse. It would be interesting to find solution of the field equation 
with radiating fluid and electromagnetic field in order to predict collective 
role of the electromagnetic field and radiation flux for the existence of the 
gravitational waves. 

The Darmois criteria for the smooth matching of the cylindrically sym- 
metric charged perfect fluid solution to an electro-vacuum charged static 
cylindrically symmetric solution leads to the following consequences: (a) 
boundary surface represents a cylinder with constant proper unit length and 
behaves as the boundary of the charged perfect fluid distributed per unit 
length of the cylinder; (b) the gravitational and Coulomb forces of the sys- 
tem balance each other on the boundary surface E. 
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